\D-A245  828 


iiRitiiiii; 


A  Second-Order  Bulk 
Boundary  Layer  Model 


David  A.  Randall  ^ 


Qingqiu  Shao 


Chin-Hoh  Moeng  2. 3 

fjo  a 

Accepted  for  publication  in 
Journal  of  the  Atmospheric  Sciences 


DTfC 

rr  i  c  T  tr 

FEB  10  1992 


Revised 

Monday,  January  6,  1992 


>  -7. -T  v  be^^n  cppiovad  I 

--1 'aUc  -  i-  _ 


^  Department  of  Atmospheric  Science 
Colorado  State  University 
Fort  Collins,  Colorado  80523 

^  National  Center  for  Atmospheric  Research 
Box  3000,  Boulder,  Colorado  80307 


^  The  National  Center  for  Atmospheric  Research  is  sponsored  by  the  National  Science  Foundation. 


ABSTRACT 


Bulk  mass  flux  models  represent  the  large  eddies  that  are  primarily  responsible  for  the 
turbulent  fluxes  in  the  planetary  boundary  layer  as  convective  circulations,  with  an  associated 
convective  mass  flux.  In  order  for  such  models  to  be  useful,  it  is  necessary  to  determine  the 
fractional  area  covered  by  rising  motion  in  the  convective  circulations.  This  fraction  can  be  used  as 
an  estimate  of  the  cloud  amount,  under  certain  conditions.  We  have  developed  "matching" 
conditions  that  relate  the  convective  mass  flux  to  the  ventilation  and  entrainment  mass  fluxes. 
These  are  based  on  conservation  equations  for  the  scalar  means  and  variances  in  the  entrainment 
and  ventilation  layers.  Methods  are  presented  to  determine  both  the  fractional  area  covered  by 
rising  motion  and  the  convective  mass  flux.  We  also  show  how  the  requirement  of  variance 
balance  can  be  used  to  relax  the  "well  mixed"  assumption.  The  vertical  structures  of  the  mean  state 
and  the  turbulent  fluxes  are  determined  analytically.  Several  aspects  of  our  simple  model's 
formulation  arc  evaluated  using  results  from  large-eddy  simulations. 
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1.  Introduction 


The  "bulk"  approach  to  parametric  representation  of  boundary-layer  processes  in  large- 
scale  models,  pioneered  by  Deardorff  (1972)  and  further  developed  by  Randall  (1976),  Benoit 
(1976),  and  Suarez  etal.  (1983)  involves  a  simple  planetary  boundary  layer  (PBL)  model  in  which 
some  aspects  of  the  vertical  structure  of  the  mean  state  are  parameterized.  Among  the  parameters 
introduced  to  represent  the  mean  state  are  the  PBL  depth,  which  is  prognostically  determined,  and 
"jumps"  or  discontinuities  at  the  PBL  top.  The  use  of  jumps  amounts  to  a  concession  that, 
although  the  fine  structure  near  the  PBL  top  is  important  for  the  PBL  physics,  it  is  unresolvable  by 
any  grid  that  can  be  used  in  a  large-scale  model.  Extensive  results  from  a  PBL  parameterization 
based  on  a  bulk  model  have  been  reponed  by  Randall  et  al.  (1985). 

The  key  shortcomings  of  existing  bulk  models  are  their  assumption  of  vertical  homogeneity 
for  conservative  variables  and  their  inability  to  predict  fractional  cloud  amounts.  In  attempts  to 
address  these  dericiencies,  a  number  of  recent  studies  (cited  below)  have  made  use  of  a  convective 
mass  flux  parameterization. 

The  mass  flux  concept  was  invented  for  use  in  a  cumulus  parameterization  by  Arakawa 
(1969).  It  was  later  adapted  to  the  problem  of  boundary  layer  parameterization  by  Betts  (1973, 
1983),  Albrecht,  et  al.  (1979),  Hanson  (1981),  Penc  and  Albrecht  (1986),  Wang  and  Albrecht 
(1986,  1990),  Randall  (1987),  and  Chatfield  and  Brost  (1987).  We  refer  to  bulk  boundary  layer 
models  that  use  the  mass  flux  concept  as  "bulk  mass  flux  models."  In  such  models,  the  convective 
mass  flux  is  assumed  to  be  associated  with  convective  circulations  that  have  ascending  and 
descending  branches.  Several  of  the  modeling  studies  just  cited  allow  the  possibility  that 
cloudiness  can  occur  (or  not)  in  either  branch.  The  concept  of  "convective  circulations"  has  also 
been  used  in  observational  studies  (based  on  conditional  sampling  and/or  joint  distribution 
functions)  by  Lcnschow  and  Stephens  (1980,  1982),  Greenhut  and  Khalsa  (1982),  Wilczak  and 
Businger  (1983),  Mahrt  and  Paumier  (1984),  Grossman  (1984),  Khalsa  and  Greenhut  (1985), 
Penc  and  Albrecht  (1986),  and  Young  (1988  a,  b).  Recently,  it  has  been  applied  to  analyze  the 
results  of  large  eddy  simulations,  by  Schmidt  and  Schumann  (1989),  Schumann  and  Moeng  (1991 
a,b)  and  Moeng  and  Schumann  (1991). 

The  existing  bulk  mass  flux  models  do  not  match  the  fluxes  associated  with  the  convective 
mass  flux  with  those  driven  by  ventilation  at  the  surface  and  entrainment  at  the  top.  Also,  no 
existing  bulk  mass  flux  model  includes  a  physically  based  method  to  determine  <t,  the  fractional 
area  covered  by  rising  motion.  Finally,  most  of  the  existing  bulk  mass  flux  models  have  retained 
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the  "well  mixed"  assumption,  with  the  notable  exceptions  of  Betts  (1973),  Albrecht  et  al.  (1979), 
and  Wang  and  Albrecht  (1990). 

This  paper  addresses  all  three  of  these  problems. 

2.  Convective  Mass  Flux  Model 

The  basic  framework  of  our  model  is  shown  in  Fig.  1.  The  level  just  above  the  PBL  top  is 
denoted  by  subscript  B+,  while  the  Eanh’s  surface  is  denoted  by  S-.  We  define  an  infinitesimal 
"ventilation  layer"  just  above  the  Earth’s  surface,  and  an  uifinitesimal  "entrainment  layer"  just 
below  the  PBL  top.  These  are  indicated  by  stippling  in  Fig.  1.  The  ventilation  layer  is  more 
conventionally  known  as  the  surface  layer.  The  entrainment  layer  is  the  region  within  which  the 
turbulent  fluxes  drop  sharply  from  finite  values  to  zero.  Caughey  et  al.  (1982)  and  Nicholls  and 
Turton  (1986)  described  the  entrainment  layer  as  a  thin  region  of  weak  organized  vertical  motions 
and  vigorous  small-scale  mixing.  The  top  of  the  ventilation  layer  will  be  denoted  by  subscript  S, 
and  the  base  of  the  entrainment  layer  by  subscript  B.  The  depth  of  the  PBL  (in  terms  of  pressure) 
is  denoted  by  dpsi- 

The  generic  variable  yf  will  be  used  to  represent  a  prognostic  intensive  scalar  such  as  the 
dry  static  energy,  the  mixing  ratio  of  water,  or  a  velocity  component.  Area-averaged  values  of  yr 
arc  denoted  by  y/.  The  upward  turbulent  flux  of  yf  is  denoted  by  Fy, .  The  mid-level  of  the  PBL  is 
indicated  by  the  dashed  line  in  Fig.  1;  it  is  representative  of  the  interior  of  the  PBL,  and  will  be 
denoted  by  subscript  "I". 

An  entrainment  mass  flux,  E,  carries  mass  across  the  PBL  top,  and  is  closely  related  to  the 
turbulent  fluxes  at  the  base  of  the  entrainment  layer.  Correspondingly,  a  ventilation  mass  flux,  V, 
is  associated  with  the  surface  fluxes;  in  conventional  parlance,  V  is  the  product  of  the  surface  air 
density,  the  surface  wind  speed,  and  a  transfer  coefficient  The  fluxes  at  the  top  of  the  ventilation 
layer  are  assumed  to  be  approximately  equal  to  those  at  the  surface.  As  explained  later,  the  model 
incorporates  diagnostic  balances  for  mass,  y/,  and  yf'^ ,  for  both  the  entrainment  and  ventilation 
layers.  These  layers  are  assumed  to  be  thin  enough  so  that  these  approximations  are  applicable. 

Within  the  ventilation  layer,  the  turbulent  fluxes  have  to  be  carried  by  small  eddies,  since 
the  organized  vertical  motions  associated  with  the  convective  circulations  must  vanish  there.  Within 
the  entrainment  layer  the  organized  vertical  motions  associated  with  the  convective  circulations 
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become  negligible.  Since  the  turbulent  fluxes  vanish  above  the  PBL  top,  however,  it  is  not 
necessary  to  hypothesize  that  small-eddy  fluxes  are  important  in  the  entrainment  layer. 
Nevertheless,  the  existence  of  small  eddies  within  the  entrainment  layer  is  well  known  from 
observations  (e.g.,  Rayment  and  Readings,  1974). 

As  indicated  in  Fig.  1,  we  assume  that  the  turbulent  fluxes  in  the  interior  of  the  PBL  are 
entirely  due  to  the  convective  circulations,  with  rising  branches  covering  fractional  area  a,  and 
sinking  branches  covering  fractional  area  I  -  a.  Some  authors  (e.g.  Greenhut  and  Khalsa,  1982) 
have  suggested  a  third,  "environmental"  domain  in  which  the  vertical  motion  is  nearly  zero.  We  are 
not  willing  to  accept  such  a  complication  without  strong  evidence  that  it  is  really  necessary.  The 
results  of  Schumann  and  Moeng  (1991  a)  suggest  that  it  is  not 

There  have  been  numerous  observations  and  numerical  simulations  yielding  values  of  a, 
based  on  various  sampling  methods.  The  methods  used  and  results  obtained  are  summarized  in 
Table  1.  Note  that  the  various  studies  are  based  on  data  from  several  different  PBL  regimes,  and 
employ  several  different  definitions  of  a.  They  have  yielded  a  variety  of  numerical  values  for  G. 

For  example,  Manton  (1977),  Coulman  (1978),  and  Lcnschow  and  Stephens  (1978)  have  adopted 
definitions  involving  not  only  the  vertical  velocity  fluctuations  but  also  the  thermodynamic 
fluctuations.  In  these  three  studies,  (Tis  found  to  be  considerably  less  than  0.5.  All  of  the 
remaining  studies  have  adopted  definitions  based  only  on  the  sign  of  the  vertical  velocity.  Greenhut 
and  Khalsa  (1982)  and  Schumann  and  Moeng  (1991  a)  defined  three  domains:  updrafts, 
downdrafts,  and  an  environment  The  updrafts  consisted  of  regions  in  which  the  vertical  velocity 
exceeded  a  positive  threshold.  If  we  take  the  fractional  area  covered  by  updrafts  to  be  cr,  then 
Greenhut  and  Khalsa  obtained  <7=  0.16,  by  far  the  smallest  value  obtained  in  any  of  the  studies, 
while  Schumann  and  Moeng  obtained  <T=  0.35.  We  note,  however,  that  Schumann  and  Moeng 
concluded  that  little  is  to  be  gained  by  defining  an  “environment”  for  the  updrafts  and  downdrafts. 

The  remaining  three  studies  listed  in  Table  1  [i.e.  those  by  Lamb  (19787,  Young  (1988  a). 

and  Nicholls  (1989)]  adopted  definitions  of  <T  that  arc  essentially  the  same  -s  ours,  namely  that  a - 

is  the  fractional  area  covered  by  positive  vertical  velocity.  The  first  two  authors  obtained  values  of  ^ 
G  equal  to  or  slightly  less  than  0.5.  Only  the  observations  of  Nicholls  (1989),  obtained  by  aircraft 
in  stratocumulus  clouds,  show  G  near  0.7.  We  conclude  that  c  is  by  no  means  an  universal 
constant,  but  that  with  the  definition  we  have  adopted  here  we  should  expect  a  to  be  not  too  : — 
different  from  1/2  under  some  conditions  of  interest 


Consider  an  arbitrary  scalar  if/,  satisfying  a  conservation  equation  of  the  form 


Statement  A  per  teiecon 
Dr.  Robert  Abb^y 
.ALrlinp.ton,  VA  22217  5 

2/7/92 


1^/ 


L(py;)  =  -V{pVy/)--^{pw\i/)  +  S^, 


where  p  is  the  density,  which  is  quasi-constant  in  time  and  the  horizontal,  as  in  the  usual  anelastic 
approximation;  V  is  d  -  horizontal  velocity  vector,  w  is  the  vertical  velocity;  and  S^f  is  the  source  of 
Vr  per  unit  mass  per  unit  time.  The  local  time  derivative  and  the  V  operator  are  defined  on  constant 
height  surfaces. 


Area  averages  satisfy 


Here  subscripts  u  and  d  denote  upward  and  downward  moving  parcels,  respectively.  The  fluxes 
associated  with  the  convective  circulations  are  given  by 

=  pw>'  p[(w„  -  vv)(  -  vr)a  +  (wj  -  w)(  y/j  -  y/){l  -  (j)] 
^Mc{Vu-¥d)> 


where 


M,  =  pa{l-a){w^-wj) 


is  the  convective  mass  flux. 


For  later  use,  we  note  that  the  variance  of  y/  is  given  by 

¥'^=(f{¥u-¥)  +U-<yi¥d-¥) 
=  a{l-<T){yf^-y/jf 


while  the  plume-scale  variance  transport  can  be  written  as 
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pw'^l/Y  =  p[^(^u-w)(¥u-¥f  +  U-<y)(^d-^)(¥d-¥) 
=  pa{I-  a){l  - 2<t)(w„  -  w^X 

x2 


=  {l-2a) 


ilA 


(2.6) 


According  to  (2.6),  pw'y/Y  is  positive  (upward  variance  transport)  if  dis  less  than  1/2, 
and  negative  (downward  variance  transport)  if  dis  greater  than  1/2.  For  d=  112,  the  variance 
transport  vanishes.  It  follows  that  if  the  ventilation  layer  is  exporting  variance  [production  exceeds 

dissipation,  and  (pn'VV')^  is  upward]  then  ds  must  be  less  than  1/2;  this  would  normally  be  the 
case  for  a  clear  convective  PBL  driven  by  surface  heating.  Similarly,  if  the  entrainment  layer  is 
exporting  variance  [production  exceeds  dissipation,  and  [pw'y/'ij/'^^  is  downward]  then  de  must 
be  greater  than  1/2. 

Now  and  firom  time  to  time  throughout  the  rest  of  this  paper,  we  investigate  the  consistency 
of  our  simple  model  with  LES  results  obtained  by  C.-H.  Moeng.  These  results  are  similar  to  the 
stratocumulus  simulations  of  Moeng  (1986),  except  that  the  domain  size  was  doubled  to  5  km  x  5 
km  X  1  km,  and  80^  grid  points  were  used.  The  vertical  resolution  was  12.5  m.  The  horizontal 
resolution  was  62.5  m.  We  have  used  six  LES  history  records,  spaced  250  simulated  seconds 
apart.  These  same  LES  results  were  analyzed  by  Schumann  and  Moeng  (1991  a;  their  "STBL " 
case). 


Fig.  2  shows  h/Cp  and  q,,  plotted  as  functions  of  height  Here  hsSy-Lqg,  Cp  is  the 

specific  heat  of  air  at  constant  pressure,  and  qi  is  the  total  mixing  ratio  (vapor  plus  liquid).  We  use 
5v  to  denote  the  virtual  dry  static  energy,  and  to  denote  the  liquid  water  mixing  ratio.  The  figure 
shows  averages  over  the  six  LES  history  records.  The  expected  well  mixed  interior  and  near¬ 
discontinuities  at  the  top  and  bottom  of  the  PBL  are  readily  apparent 

We  propose  a  new  method  to  determine  (land  Me,  based  entirely  on  the  vertical  velocity 
statistics,  and  suggested  by  the  analysis  of  Moeng  and  Rotunno  (1990).  The  idea  is  to  use  (2.2), 
(2.5),  and  (2.6)  with  =  w,  and  to  solve  these  three  equations  for  the  three  unknowns  <t,  and 
wj.  The  convective  mass  flux  can  then  be  evaluated  using  (2.4).  We  find  that 
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where 


(2.7) 


a 


2a^4  +  S, 


2’ 


(2.8) 

is  the  skewness  of  the  vertical  velocity.  According  to  (2.7),  or  is  less  than  1/2  when  the  skewness 
is  positive,  and  greater  than  1/2  when  the  skewness  is  negative.  The  form  of  (2.7)  guarantees  that 
0  <  O'  <1.  The  convective  mass  flux  satisfies 


(2.9) 

One  advantage  of  this  approach  is  that  it  guarantees  consistency  among  the  mass  flux  and  the 
variance  of  the  vertical  velocity,  within  the  context  of  the  mass  flux  model.  A  second  advantage  is 
that  it  can  be  used  to  determine  crand  Me  a  higher-order  closure  model.  Of  course,  we  cannot 
use  this  method  to  determine  o  and  in  a  bulk  mass  flux  model,  because  such  a  model  does  not 
provide  the  needed  input  data,  e.g.  is  not  available.  Later,  in  Section  4,  we  present  a  method 
that  can  be  used  to  determine  a  and  in  a  bulk  mass  flux  model. 

Fig.  3  a  shows  the  profiles  of  and  obtained  by  the  method  described  above,  and 
also  by  directly  sampling  the  LES  results.  The  profiles  are  quite  similar,  but  the  method  proposed 
here  systematically  overestimates  the  magnitudes  of  and  wd .  Fig.  3  b  shows  the  profiles  of  <J 
and  Me,  obtained  from  the  LES  results  by  the  method  described  in  the  preceding  paragraph. 
Although  a  is  close  to  0.5  at  all  levels,  it  decreases  noticeably  upward  through  most  of  the  PBL's 
depth,  from  about  0.6  near  the  surface  to  about  0.4  slightly  below  the  PBL  top.  This  is  consistent 
with  the  analysis  of  Moeng  and  Rotunno  (1990).  The  convective  mass  flux  has  a  roughly  parabolic 
profile,  with  a  maximum  slightly  above  the  mid-level  of  the  PBL,  and  much  smaller  values  near  the 
surface  and  the  PBL  top.  The  maximum  value  of  Me  is  about  0.45  kg  m*2  s'L 

Using  (2.2)  and  (2.3),  we  can  evaluate  yfu  and  y/d  as  follows: 
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—  /-<T  „ 

Me 

(2.10) 

Me 

(2.11) 

The  values  of  \ffu.  and  Wd  obtained  directly  by  conditional  sampling  (based  on  the  sign  of  the 


vertical  velocity)  are  slightly  different  from  those  obtained  from  the  fluxes  using  (2.10  -  1 1).  This 
is  shown  in  Fig.  4  a  for  htCp  and  in  Fig.  4  b  for  q, .  In  these  figures,  the  solid  lines  show  the 
updraft  and  downdraft  properties  obtained  from  (2.10  -  11),  and  the  dashed  lines  show  the 
corresponding  values  based  on  conditional  sampling.  The  updrafts  are  "warm"  and  "wet,"  and  the 
downdrafts  are  "cold"  and  "dry."  To  give  the  correct  fluxes,  (2. 10  - 11)  make  the  updrafts  warmer 
and  wetter,  and  the  downdrafts  colder  and  drier,  than  the  corresponding  updraft  and  downdraft 
properties  obtained  by  conditional  sampling  (c.f.  Young,  1988  a;  Wang  and  Albrecht,  1990; 
Schumann  and  Moeng,  1991  a).  The  discrepancies  are  fairly  small,  however,  compared  to  the 
differences  between  the  updraft  and  downdraft  properties.  These  discrepancies  are  due  to  the  "top 
hat"  profiles  assumed  in  the  convective  mass  flux  model.  In  a  sense,  the  mass  circulation  of  the 


convective  mass  flux  model  produces  fluxes  less  efficiently  than  the  more  realistic  simulation 
produced  through  LES.  In  the  remainder  of  this  paper,  we  use  \ffu  and  xffd  as  determined  from  (2. 10 


-  11). 


3.  Matching  the  fluxes 


The  surface  fluxes  are  assumed  to  satisfy  the  usual  bulk  aerodynamic  formula, 

(3.1) 

where  V  is  the  "ventilation  mass  flux,"  which  is  usually  written  as  the  product  of  surface  wind 
speed,  a  transfer  coefficient,  and  the  surface  air  density. 
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We  assume  that  the  fluxes  at  the  top  of  the  ventilation  layer  are  entirely  due  to  the 
convective  circulations,  and  that  the  small-eddy  fluxes  are  negligible  there.  Using  this  assumption 
with  (2.3),  (3.1),  and  also  using  our  assumption  that  the  ventilation  layer  is  thin,  we  can  write 

V(¥s--¥s)  =  ^c.s(¥u'¥d)s- 

(3.2) 

This  condition  implies  a  consistency  between  the  fluxes  obtained  from  the  bulk  formula  and  those 
determined  from  the  mass  flux  model.  Wang  and  Albrecht  (1990)  did  not  impose  (3.2)  or  the 
corresponding  condition  at  the  PBL  top  (discussed  later). 

At  level  S,  the  parcels  rising  away  from  the  lower  boundary  must  be  "charged"  with  the 
properties  of  the  boundary.  We  cannot  assume,  however,  that  the  properties  of  the  updrafts  at  level 
S  are  the  same  as  those  of  the  boundary,  because  there  can  be  very  strong  gradients  across  the 
ventilation  layer.  The  small  eddies  of  the  ventilation  layer  rapidly  dilute  air  that  has  been  in  conuct 
with  the  boundary,  by  mixing  it  with  air  that  has  recently  descended  from  the  interior  of  the  PBL. 

As  a  result,  ¥d)\>  from  (3.2),  this  implies  that  V.  In  order  to  take 

this  into  account,  we  introduce  a  nondimcnsional  parameter,  Xv,  such  that 

(¥u)s-¥s  =  Xv(¥s--¥sy^ 

(3.3) 

in  case  =  /,  we  get  (iffJs  =  ¥s-  •  Smaller  values  of  indicate  stronger  mixing  by  the  small 
eddies  of  the  ventilation  layer.  We  expect  0  <  Xv  «  combining  (2.2),  (3.2),  and  (3.3),  we 


find  that 


;fv^  =  l-(Ts, 


(3.4) 


where  Gy  =  5  /  V  .This  is  a  "continuity  equation"  for  the  eddies,  expressing  a  relationship 

between  the  convective  mass  flux  and  the  ventilation  mass  flux. 


In  the  preceding  discussion,  it  has  been  tacitly  assumed  that  Xv  is  independent  of  i.  e.. 
that  a  single  "mixing"  parameter satisfies  (3.3)  whether  y^is  moist  static  energy,  total  water 
mixing  ratio,  or  some  other  intensive  scalar.  This  assumption  is  supported  by  (3.4);  if  1^.  and 
Os  are  all  independent  of  ¥>  Xv  must  also  be  independent  of  y/.  Such  independence  suggests 
that  Xv  is  a  useful  concept 
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We  now  apply  a  similar  analysis  to  the  entrainment  layer.  The  assumption  that  it  is  thin 
yields  the  familiar  "jump"  relation  between  and  tl.e  entrainment  rate: 

-a* 

{F^)b  =  -  v^a)  -  dz. 

•• 

(3.5) 

Here  E  is  the  rate  at  which  mass  is  entrained  across  the  PBL  top.  In  (3.5),  we  follow  Lilly  (1968) 
by  keeping  the  term,  which  represents  a  possible  concentrated  entrainment-layer  "source"  of  y/, 
(e.g.,  due  to  radiation).  We  now  assume  that  the  fluxes  at  the  base  of  the  entrainment  layer  are 
entirely  due  to  the  convective  circulations,  and  that  the  "small-eddy"  fluxes  are  negligible.  Then,  by 
comparing  (2.3)  and  (3.5),  and  using  (2.2),  we  obtain 


Mc,B(¥d-¥)B  =  E:(yB(¥B+-  ¥b)  +  <^b 

-’fl 

(3:6) 

At  this  point,  we  introduce  a  mixing  parameter  Xe^  by  analogy  with  (3.3).  To  allow  for  the 
effects  of  the  concentrated  source,  however,  we  include  an  additional  term: 

(¥,-¥^b^Xe(¥b^-¥b)-^^  jKdz- 

(3.7) 

Here  A  is  a  coefficient  that  is  determined  below.  According  to  (3.7),  the  properties  of  the 
descending  air  at  level  B  are  related  to  those  of  the  free  atmosphere  just  above  the  PBL  top,  as 
modified  by  small-eddy  mixing  and  the  effects  of  any  concentrated  source  within  the  entrainment 
zone.  Since  there  is  a  sharp  gradient  of  across  the  entrainment  layer,  we  expect  0<  Xe«  !■ 

The  mixing  parameter  ;|f£is  closely  related  to  the  parameter  x  discussed  by  Siems  et  al. 
[1989;  see  also  Albrecht  etal.  (1985)  and  NichoUs  and  Turton  (1986)].  We  can  interpret  Xe  as  the 
value  of  X  associated  with  the  downdraft  air  at  level  B.  Further  discussion  is  given  later  in  this 
Section. 

Comparing  (3.6)  with  (3.7),  we  find  that 
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(3.8) 


(-^c.bXe-^  B*-  +  +  (ya)  = 

-s 


In  case  the  source  term  of  (3.8)  vanishes,  we  obtain 


XeG^  -  <Jb, 

(3.9) 

where  =  M^  g  /  E .  This  relationship  does  not  involve  y/;  it  must,  therefore,' apply  for  all  y/.  To 
ensure  that  (3.9)  will  be  satisfied  even  when  the  source  term  of  (3.8)  is  not  zero,  we  must  choose 

X  =ag/ Mc.b- 

(3.10) 

We  can  interpret  (3.9)  as  another  "continuity  equation",  analogous  to  (3.4).  Again,  we 
have  tacitly  assumed  that  ^ffis  independent  of  the  species  under  consideration.  This  assumption  is 
consistent  with  (3.9),  since  Mc,«,  E,  and  Gg  are  independent  of  species.  We  can  use  (3.9)  to 
eliminate  Me,  a  in  (3.6),  or,  alternatively,  use  (3.10)  to  eliminate  X  in  (3.7);  either  way,  the  result 
is: 


{¥d)B  =  XE¥B*-^(^’'XE)¥B-^^  ]Xdz. 

^  -8 

(5.11) 

According  to  (3.1 1),  the  descending  air  at  level  B  has  the  properties  of  the  free  atmosphere,  except 
as  modified  by  mixing  (when  Xe  <  O  and  by  the  concentrated  source.  Caugbey  et  al.  (1982)  and 
Nicholls  (1989)  have  reported  observations  of  cool  downdrafts  in  the  upper  portions  of 
stratocumulus  cloud  sheets.  They  concluded  that  the  sinking  air  had  been  radiatively  cooled  near 
the  cloud  top.  Such  effects  are  represented  by  the  S^term  of  (3.1 1). 

This  term  is  inversely  proportional  to  E,  which  means  that  radiative  cooling  in  the 
entrainment  layer  can  produce  negatively  buoyant  parcels  most  effectively  if  the  entrainment  rate  is 
slow.  This  suggests  that  entrainment  driven  by  radiative  cooling  near  cloud-top  tends  to  be  self- 
limiting. 

According  to  (3.4),  ventilation-layer  dilution  becomes  more  effective  (in  other  words.  Xv 
decreases)  as  the  convectiv'*  mass  flux  increases  relative  lo  V  ( I  -  Gs  ).  The  ventilation  mass  flux 
times  the  fractional  area  covered  by  the  incoming  downdrafts  is  a  measure  of  the  rate  at  which  the 
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updrafts  leaving  the  ventilation  layer  can  be  supplied  with  air  that  has  been  charged  with  surface 
properties,  and  the  convective  mass  flux  is  a  measure  of  the  rate  at  which  this  air  is  removed  from 
the  surface  layer.  The  stronger  the  convective  mass  flux  becomes,  the  less  effectively  ventilation- 
layer  air  can  be  charged  with  surface  properties  before  it  is  carried  away  into  the  interior  of  the 
PBL.  A  similar  interpretation  can  be  given  for  (3.9). 

If  we  combine  (3.4)  and  (3.9),  and  assume  that  a  is  independent  of  height  through  the 
depth  of  the  PBL,  so  that  <75  =  (Tg,  we  find  that 


a 


1  + 


_l _ 

Xe^e 


(3.12) 

Suppose  that  Me  is  also  independent  of  height,  so  that  Mc,s  ~  Mcjb-  If  we  substitute  (3.12)  back 
into  either  (3.4)  or  (3.9),  and  use  the  definitions  of  Gy  and  Gf,  we  obtain 


w  _  {EIxeIVIXv) 

'  {^ixEhiyixvY 


With  the  assumption  that  Me  is  independent  of  height,  (3.12)  can  be  rewritten  as 


(3.13) 


<7  = 


(3.14) 

Of  course,  bulk  PBL  models  normally  include  parameterizations  for  V  and  E.  If  we  have  a 
parameterization  for  then  (3.14)  can  be  used  to  determine  a.  Note  that  the  forms  of 

(3.12)  and  (3.14)  ensure  thatO  <  a<  I.  We  use  both  (3.13)  and  (3.14)  later  in  this  paper. 


The  profiles  of  crand  Me  shown  in  Fig.  3  are  not  independent  of  height.  Nevertheless,  the 
assumption  that  these  two  parameters  are  vertically  uniform  is  of  interest  as  a  particularly  simple 
special  case. 


Now  we  test  the  main  conclusions  of  this  Section,  using  the  LES  results. 
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First  we  have  to  identify  levels  S-,  S,  B  and  B+  based  on  the  LES  output.  In  this  we  are 
guided  not  only  by  the  mean  structures  but  also  by  the  turbulence  results,  since  we  have  defined 
the  ventilation  and  entrainment  layers  in  terms  of  the  relative  importance  of  small -eddy  fluxes.  Fig. 
5  shows  the  mean  structure  of  the  total  moisture  qt  (g  kg'l.  solid  line)  and  the  ratio  of  the  subgrid 
moisture  flux  to  the  total  moisture  flux.  The  ratio  becomes  appreciable  only  near  the  top  and  bottom 
of  the  PBL.  We  identify  the  entrainment  and  ventilation  layers  as  those  regions  within  which  the 
ratio  exceeds  0.1.  We  have  also  taken  into  account  the  mean  structure.  The  entrainment  and 
ventilation  layers  that  we  have  identified  are  indicated  in  Fig.5.  The  height  of  level  S-  is  zero  (the 
surface),  while  level  S  is  at  62.5  m,  level  B  is  at  475  m,  and  level  B-i-  is  at  537.5  m.  The  average 
cloud  top  is  near  504  m.  Since  the  vertical  grid  spacing  of  the  LES  model  is  12.5  m,  there  are  four 
model  layers  inside  the  ventilation  layer,  and  four  more  inside  the  entrainment  layer. 

Averaging  over  six  LES  history  records,  we  obtain  the  following  numerical  values  for  the 
mean-state  variables: 


295.73  K, 

285.45  K. 

4.22gkg~[ 

II 

7.18gkg-{ 

285.51  K, 

^S-  ^^p~ 

288.67  K, 

11 

7.23gkg-{ 

9.80  gkg~y 

These  values  are  indicated  in  Fig.  2. 

We  have  evaluated  the  ventilation  mass  flux  from  the  LES  results,  as  the  ratio  of  the  surface 
flux  of  h  or  qt  to  the  surface-air  difference  in  the  same  quantity,  as  indicated  in  (3.1).  Averaging 
over  six  LES  history  records,  the  ventilation  mass  flux  based  on  h  is  1.04  x  10'^  kg  m‘2  s■^  while 
that  based  on  qt  is  0.965  x  10*2  kg  m*2  s**.  We  estimate,  therefore,  that  the  correct  value  of  V  is 
about  1.0  X  10*2  kg  m*2  s*L 
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Table  2  gives  the  entrainment  mass  flux,  E,  as  evaluated  from  the  LES  data  by  the 
following  methods: 

Method  HA: 


Method  HB: 


*  Ah 


Method  QA: 


Method  QB: 


Method  Z: 


Ah 


-pWB 


Here  zt  is  the  cloud  top  height  Methods  HA  and  QA  include  the  small  contributions  from  the 
fluxes  at  level  B+,  and  so  are  a  priori  more  reliable  than  methods  HB  and  QB.  Our  results  indicate 
that  the  effect  of  fluxes  at  level  B+  on  the  estimated  entrainment  rate  amounts  to  about  20%. 
Method  Z  does  not  give  a  reliable  value  for  the  entrainment  rate  because  of  uncertainties  in 
determining  the  values  of  27  and  especially  its  time  change;  we  have  shown  only  the  average  value 
based  on  the  total  length  of  the  record  available  to  us.  From  the  results  given  in  Table  2.  we 
conclude  that  the  correct  value  of  E  is  about  0.82  x  lO'^  kg  m'2  s’T  This  corresponds  to  an 
“entrainment  velocity”  of  about  25  m  hr*  or  about  600  m  day*,  which  is  in  line  with  observations 
(e.g.,  Nicholls  and  Turton,  1986).  Note  that  E  and  V  are  quite  comparable  in  magnitude. 

We  have  argued  that  Xv  and  xe  are  independent  of  species.  Do  the  LES  results  support 
this?  Substituting  (2.10)  and  (2.11)  into  (3.3)  and  (3.6),  and  using  (3.10),  we  obtain  expressions 
for  XV  Sind  XE  - 
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Xy  = 


1 - 

1 

\ _ ' 

^  ’^cS  j 

¥s--¥s 

(3.15) 


(3.16) 


For  y/  =  qi,  we  have  S^=0,  while  for  yr  =  we  have 


(3.17) 


where 


^  Jrad 


is  the  heating  rate  due  to  longwave  radiation. 


Using  the  methods  of  Section  2,  we  have  evaluated  a  and  the  convective  mass  fluxes  at 
levels  B  and  S.  The  results,  averaged  over  the  six  history  records,  are  as  =  0.595,  <Jb  =  0.352. 
Mc,s  =  0.316  kg  m*2  s‘^  and  Mc,b  -  0.235  kg  m*2  s*^  From  the  mass  fluxes,  we  obtain  Gy  = 
31.56,  and  Ge  =  28.66.  We  have  also  evaluated  AR,  by  using  (3.17).  The  result  is  zW?  =  82.1  W 
m'2.  In  a  similar  way,  we  have  determined  that  the  radiative  flux  "jump"  across  the  ventilation  layer 
is  a  warming,  of  2.2  W  m*2. 

We  have  also  determined  Xv  and  xe  using  each  of  the  six  LES  history  records.  For  this 
purpose,  we  used  (3.3)  and  (3.7),  with  the  numerical  values  of  y/d,  and  y?  given  earlier.  The 
results  are  given  in  Table  3.  The  XE  s  obtained  for  yr  =  /i  and  y/=  qt  agree  very  well.  Note  that 
Xy.h  is  generally  a  little  bit  larger  than  xv.^  If  we  include  the  small  radiative  flux  jump  across  the 
ventilation  layer,  we  find  that  the  two  xv 's  agree  as  well  as  the  two  XE  s  do.  On  the  basis  of  these 
results,  we  conclude  that  xv=  1-22  x  10*2,  and  that  Xe  =  102  x  10*2.  jhe  good  agreement 
between  the  values  of  x  obtained  with  h  and  qt  suppons  our  assertion  that  Xv  ^ud  xe  are 
independent  of  species.  As  expected,  both  xv  and  xe  arc  small  compared  to  one.  Their  numerical 
values  are  quite  similar. 

Since  we  have  chosen  V  and  E  so  that  (3.1)  and  (3.5)  are  satisfied,  the  good  results 
presented  above  for  Xv  and  Xe  ensure  that  (3.4)  and  (3.9)  will  be  satisfied.  Substituting  our  "best 
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estimates”  of  E.  V,  xv  and  into  (3.14),  we  find  that  the  effective  height-independent  value  of  a 
is  0.505.  in  reasonably  good  agreement  with  the  LES  results  shown  in  Fig.  3. 

Using  the  numerical  values  given  above  for  E,  V,  xe>  Xv,  we  find  from  (3.13)  that  the 
effective  height-independent  value  of  Me  is  0.406  kg  m'2  s'^.  This  is  slightly  larger  than  the  value 
of  MeS  obtained  directly  from  the  LES  results,  and  is  substantially  larger  than  the  value  of  Mc^. 
It  is  fairly  close  to  the  maximum  value  of  Me  shown  in  Fig.  3. 


Before  leaving  this  Section,  we  offer  an  alternative  interpretation  of  xe  ■  The  air  descending 
at  level  B  is  a  mixture  of  updraft  air  that  has  passed  through  the  entrainment  layer,  and  newly 
entrained  air  from  level  B+.  Let/  be  the  fraction  of  air  from  level  B-t-,  i.e. 


(3.18) 

For  simplicity  we  have  assumed  no  concentrated  sources  or  sinks;  note,  however,  that  this 
assumption  does  not  alter  the  result  below.  Comparing  (3.18)  with  (3.7),  we  find  that 


/  = 


Xe 

XE-^<^B(i-XE)' 


(3.19) 

Substituting  the  numerical  values  of  crand  Xe  from  above,  we  obtain  /  =  0.028.  This  means  that, 
for  the  case  studied  here,  about  3  %  of  the  air  descending  at  level  B  has  just  been  entrained:  the 
remainder  is  “recycled”  updraft  air. 

Appendix  A  presents  an  analysis  of  the  budgets  of  for  the  entrainment  and  ventilation 


layers. 


4.  The  convective  mass  flux  and  the  turbulence  kinetic  energy 


In  this  Section,  we  present  methods  to  determine  crand  Me  inside  the  convective  mass  flux 
model,  using  the  turbulence  kinetic  energy  (TKE)  and  the  entrainment  and  ventilation  mass  fluxes 
as  inputs. 

Obviously,  there  has  to  be  a  close  connection  between  the  convective  mass  flux  and  the 
turbulence  kinetic  energy  (TKE).  Let  aj  be  the  fraction  of  the  vertically  integrated  TKE  that 
resides  in  the  vertical  component  of  the  motion.  We  can  write 
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dz  =  4l:*'pw'^  dz . 


Here  ^.v/  is  the  vertically  averaged  TKE  per  unit  mass.  If  the  TKE  is  equipartitioned  among  the 
three  convective  velocity  components,  then  dj  =  ll3.  Using  the  LES  results,  and  considering  only 
the  resolved-scale  motions  in  the  region  25.  <  z  ^  ,  we  find  that  e,v/  =  0.7 12  m^  s*^,  and  dj  = 

0.281. 

Using  (2.3),  we  can  express  in  terms  of  the  difference  in  vertical  velocity  between  the 
updrafts  and  the  downdrafts.  Next,  (2.4)  can  be  used  to  rewrite  this  difference  in  terms  of  the 
convective  mass  flux,  giving 

— 

w  =—5 - - - . 

p^ad-G) 

(4.2) 

Because  of  the  way  we  have  determined  a  and  Me  from  the  LES  results  (sec  Section  2),  (4.2)  has 
to  be  exactly  satisfied  when  is  evaluated  ft-om  the  LES.  By  substitution  from  (4.2),  (4.1)  can 
be  rewritten  as 


^s^MWydz  =  j/;*; 


pad -a) 


If  a  and  Me  are  independent  of  height,  then  (4.3)  reduces  to 


2aj<T(  I  - 

£'4 


The  integrals  in  (4.4)  are  easily  evaluated  for  a  given  sounding.  If  we  consider  the  density  of  the 
air  to  be  approximately  constant  with  height  in  the  PBL,  then  (4.4)  can  be  simplified  to 

M^  =  Pi^^2d^a{\-a)ef^. 


Appendix  B  gives  a  method  to  predict  cm  in  a  i^ulk  boundary  layer  model.  Then  (4.5)  can  be  used 
to  evaluate  Me  ,  provided  that  a  is  known. 
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From  the  LES  results,  we  have  already  found,  in  Section  3,  that  the  effective  heighi- 
’ndependent  values  of  Me  and  a  are  0.406  kg  s’l  and  0.505,  respectively.  Substituting 
numerical  values  into  the  rhs  of  (4.5),  we  obtain  Me  =  0.383  kg  m*2  s'^.  This  is  in  fair  agreement 
with  the  value  just  mentioned,  and  seems  reasonable  in  view  of  the  LES  results  shown  in  Fig.  3. 


By  equating  (4.5)  and  (3.13),  and  using  (3.14),  we  can  derive  a  constraint  on  the  product 
XvXe  ■  i  e. 


XvXe  = 


EV 


(4.6) 

Substituting  numerical  values  on  the  rhs  of  (4.6)  gives  XvXe  ~  L40  x  10“*.  which  should  be 
edmpared  with  the  value  inferred  in  Section  3,  i.e.  XvXe  -  (1-22  x  10'2)  x  (1.02  x  10*2)  =  1.24  x 
10-4. 


Alternatively,  we  can  set  the  Ihs  of  (4.5)  to  0.406  kg  m*2  s**,  which  is  the  value  obtained  in 
Section  3,  and  solve  (4.5)  for  aj.  This  gives  aj  =  0.316,  which  is  slightly  larger  than  the  value 
obtained  directly  from  the  LES  results,  mentioned  above.  This  small  discrepancy  in  aj  comes  from 
our  assumption  that  both  a  and  Me  arc  independent  of  height.  If  Oj  =  0.316  is  used  in  (4.6),  we 
are  guaranteed  to  obtain  the  same  value  of  XvXe  reported  in  Section  3.  Of  course,  given  a 

method  to  determine  XE  and  xv ,  we  could  use  (4.6)  to  solve  for  Oj. 

Our  plan  is  to  parameterize  XE  and  xv-  The  parameterization  must  be  consistent  with  (4.6). 
Note  that  the  denominator  on  the  rhs  of  (4.6)  contains  only  quantities  that  characterize  the  PBL  as  a 
whole,  rather  than  the  entrainment  or  ventilation  layers.  We  expect  E  and  V  to  appear 
symmetrically  in  the  expressions  for  XE  and  xv  ■  The  form  of  (4.6)  thus  suggests  the  following 
two  alternative  possibilities  for  and;KV  : 

V  E 

(4.7  a) 

or 

I  EV 

Xv-Xe-  ■>  r; — 72 —  • 

(4.7  b) 
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Obviously,  neither  (4.7  a)  nor  (4.7  b)  follows  rigorously  from  (4.6);  these  are  just  two  panicularly 
plausible  possibilities.  A  possible  motivation  for  (4.7  a)  is  that  it  pairs  V  with  xv  and  E  with  xe  • 
thus  keeping  ventilation-layer  quantities  together  and  entrainment  layer  quantities  together.  A 
motivation  for  (4.7  b)  is  that  the  product  EV  is.  in  a  sense,  characteristic  of  the  PBL  as  a  whole, 
so  that  it  can  plausibly  appear  symmetrically  in  the  expressions  for  Xe  and  Xv  ■ 

Substituting  the  previously  mentioned  numerical  values  for  the  quantities  on  the  right-hand- 
sides  of  (4.7  a),  we  obtain  Xv  -  1-31  x  lO'^  and  Xe  ~  ^  Similarly.  (4.7  b)  gives  Xv  - 

Xe  -  118  X  10‘2.  These  values  are  reasonably  close  to  those  deduced  from  the  LES  results  in 
Section  2.  i.e.  Xv  =  1-22  x  10*2  -  i.02  x  10*2. 


We  can  now  substitute  (4.7  a)  or  (4.7  b)  back  into  (3.14).  to  obtain  an  expression  for  a. 
From  (4.7  a)  we  obtain: 


a  =  f- 


(4.8  a) 


and  from  (4.7  b): 


V 


(4.8  b) 

Since  £  and  V  are  practically  equal  in  the  LES  results  under  discussion  here,  (4.8  a)  and  (4.8  b) 
give  essentially  the  same  numerical  value  for  <T,  and  both  are  reasonably  consistent  with  the  LES 
results  presented  in  Section  3,  which  gave  <7  close  to  1/2 .  It  would  be  useful  to  have  observations 
or  LES  results  for  which  £  and  V  differed  appreciably;  if  ff  was  also  known,  this  would  allow  us 
to  choose  between  (4.8  a)  and  (4.8  b). 


5.  Interiors 


Up  to  this  point  we  have  discussed  only  the  ventilation  and  entrainment  layers.  Now  it  is 
time  to  consider  the  interior  of  the  PBL.  The  purpose  of  this  Section  is  to  show  how  the  “well 
mixed”  assumption  can  be  relaxed  by  using  a  simple  second-order  closure,  within  the  context  of 
our  bulk  mass  flux  model. 


We  assume  that  the  variance  budget  of  the  PBL’s  interior  satisfies 
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p  dz 


—■^ipw'\i/'y/')-2e 
p  dz 


w 


(5.1) 

Advection  by  the  mean  flow  has  been  neglected,  and  we  have  assumed  for  simplicity  that  is  a 
conservative  variable.  (This  assumption  is  not  necessary  and  can  easily  be  relaxed.)  We  now 
model  each  term  of  (5.1),  following  the  methods  introduced  earlier.  The  variance  itself  is  replaced 
by  (2.5).  The  triple  correlation  in  the  transport  term  is  replaced  by  (2.6).  The  dissipation  is 
modeled  by 


(5.2) 

where  Tdis  is  a  dissipation  time  scale  whose  functional  dependence  on  a  will  be  discussed  later. 
With  the  assumptions  that  (Tand  Me  are  independent  of  height,  we  find  that 


d_ 

dt 


M. 


c  ) 


=  -2 


d\ff 

M,{l-2o)  1  d 

p  dz 

0(1 -a)  p  dz 

2a(l-a) 


^dis 


¥ 


(5.3) 


It  is  interesting  to  consider  the  equilibrium  (d  I  dt  =  0)  solutions  of  (5.3)  for  two  limiting 
cases.  First,  suppose  that  a  -  112  ,  so  that  the  middle  term  (the  transport  term)  of  (5.3)  drops  out. 
Then  we  get 


F  - 

^  0(1  -  a)  p  dz  ' 

This  is  a  down-gradient  diffusion  formula. 


(5.4) 


Next,  suppose  that  <t  «  7,  as  in  the  case  of  cumulus  convection.  Then  (5.3)  reduces  to 

*  ■ 

(5.5) 
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This  is  the  “compensating  subsidence”  formula  that  has  become  familiar  in  cumulus 
parameterization  theories.  (The  “detrainment”  terms  that  appear  in  cumulus  parameterizations  can 
be  included  here  too  by  relaxing  our  simplifying  assumption  that  Me  is  independent  of  height.) 

In  shon,  (5.3)  includes  both  downgradient  mixing  and  “compensating  subsidence”  as 
special  cases. 

Fig.  6  shows  the  dissipation  time  scale  obtained  from  the  LES  results,  as  a  function  of 
height,  determined  by  using 


(5.6) 

for  both  h  (solid  line)  and  qt  (dashed  line).  Near  the  mid-level  of  the  PBL,  is  on  the  order  of  a 
thousand  seconds.  It  is  considerably  shorter  near  the  surface  and  the  PBL  top.  This  might  be 
because  the  resolved-scale  variances  are  small  near  the  edges  of  the  PBL,  or  it  might  be  because 
dissipation  is  actually  more  effective  in  those  regions  where  small  eddies  are  dominant.  We  use  a 
single,  height-independent  value  of  in  our  model.  Further  discussion  of  is  given  later. 


We  can  use  these  results  in  either  of  two  ways,  both  of  which  allow  us  to  relax  the  “well 
mixed”  assumption.  On  the  one  hand,  the  equilibrium  solution  of  (5.3)  gives  a  constraint  on 
/  dz  similar  to  that  proposed  by  Wyngaard  and  Brost  (1984).  In  fact,  one  way  to  proceed  from 

here  would  be  to  adopt  the  formulae  for  and  proposed  by  Wyngaard  and  Brost  (e.g., 

their  Eq.  46),  and  to  force  consistency  with  our  model.  This  approach  would  yield  constraints  on  a 
and  Me . 


The  other  possible  application  of  (5.3)  is  as  follows.  Note  that  (5.3)  is  hyperbolic;  only  the 
first  derivative  with  respect  to  z  of  <j(I-a^F^I  M^f  appears,  and  it  is  multiplied  by  the 

“velocity”  M^l -2a)/[pa(l -<t)].  This  “velocity”  is  upward  (positive)  if  <y  <  1/2,  and 
downward  (negative)  if  <r  >  1/2.  It  follows  that  the  single  boundary  condition  on 

a(l-(7){^F^I  M^  must  be  applied  at  level  S  if  <t<  1/2,  and  at  level  B  if  (T>  1/2.  We  cannot  force 

the  surface  flux  to  satisfy  the  bulk  aerodynamic  formula  and  the  flux  at  level  B  to  satisfy  (3.5) 
simultaneously,  unless  some  additional  freedom  is  introduced. 
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The  needed  freedom  lies  in  the  choice  of  <9  i/r  /  ^,  which  is  unknown  at  this  point  anyway. 
We  now  assume  that  !  dz  is  constant  with  height,  and  drop  the  time-derivative  term  of  (5.4). 
This  leads  to  an  ordinary  differential  equation  for  f  y, ,  which,  using  hydrostatics,  can  be  written  as 


where 


{l-2a)( 

Me  {  Bp 


I  _  By/ 
j  Bp  ' 


5p*  = 


a{l-a){l-2a) ' 


(5.7) 


(5:9) 

We  assume  that  x  is  independent  of  a.  The  motivation  for  this  assumption  is  explained  in 
Appendix  C. 


The  solution  of  (5.7)  is 


{l-2a)  Bp  \  Sp* 


(5.10) 

where  A  is  an  arbitrary  constant.  In  order  to  maintain  consistency  of  (5.10)  with  the  imposed 
fluxes  at  the  top  and  bottom  of  the  PBL,  we  must  choose  appropriate  values  of  Bxj/ 1  Bp  and  A. 

In  this  way,  we  find  that 


and 


Byf  _  {l-2a) 

Bp  MeSp* 

/  -exp 

l^J 

(5.11) 
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(5.12) 

These  results  show  how  the  mean  gradient  and  the  flux  profile  are  determined  for  given  values  of 
the  fluxes  at  levels  S  and  B  and  the  other  parameters.  Note  from  (5.1 1)  that  d\t/ 1  dp  vanishes  for 
(T  =  1/2.  Inspection  of  (5.12)  confirms  that  satisfies  the  appropriate  boundary  conditions  at 
levels  S  and  B. 

Obviously,  Sp;^  /  dp»  is  a.  key  parameter  in  (5.1 1  -  12).  It  can  be  written  as 

Sp, 

(5.13) 

As  an  example.  Fig.  7  shows  how  dpf^  /  dp*  varies  with  a,  for  the  particular  case  Me  =  0.4  kg  m*2 
s*^  Spif^  =  50  mb.  The  solid  line  is  for  f  =  1  second,  and  the  dashed  line  is  for  f=  10  seconds 
(see  Appendix  C  for  a  discussion  of  t).  Note  that  dpi^  /  Sp,  passes  through  zero  for  <j  =  1/2.  As 
is  apparent  from  (5.13),  dp^^l  Sp*  decreases  as  f  increases.  In  this  example,  for  f  =  1  second 
dpi^  /  dp*  can  be  of  order  ±10^  for  a  moderately  different  from  1/2.  Smaller  values  of  ^  v/ ! 
are  favored  by  a  stronger  convective  mass  flux  and  larger  values  of  f . 

Extending  this  example.  Fig.  8  shows  how  d^f  !  dp  varies  with  a.  Again,  the  solid  line  is 
for  f=  1  second,  and  the  dashed  line  is  for  t=  10  seconds.  Here  we  consider  potential 
temperature,  6  (Fig.  8  a),  and  the  mixing  ratio  of  water  vapor,  q  (Fig.  8  b).  The  “edge  ”  fluxes 

are  assumed  to  be  =  0.1  kg  K  m'2  s**,  (Fq)^  =  -0.02  kg  K  m'^  s'*,  and  (f^)^  =  3  x  10'5  kg 

m*2  s'l  and  (f^)^  =  2.5  x  10*5  kg  m*2  s*^.  For  potential  temperature,  the  results  are  plotted  as  K 

(50  mb)*l,  and  for  moisture  they  are  plotted  as  g  kg**  (50  mb)**.  The  figure  shows  that  the  sign  of 
dd  I  ^  depends  on  <7 ;  positive  values  (upward  decrease)  occur  for  <t>  1/2,  and  negative  values 
(upward  increase)  for  <J<  1/2.  In  contrast,  the  sign  of  dq  I  dp  is  positive  (upward  decrease)  for  all 
values  of  a .  The  different  behaviors  of  dB  t  dp  and  dq  /  dp  arise  from  the  differences  in  their 
respective  presecribed  fluxes  at  levels  S  and  B.  With  the  values  of  these  fluxes  that  we  have 
prescribed  for  this  example,  observations  suggest  (e.g.,  Wyngaard  and  Brost,  1984)  that  potential 
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temperature  should  increase  upward  and  water  vapor  mixing  ratio  should  decrease  upward.  Our 
model  predicts  such  a  result  for  a  less  than  1/2.  The  numerical  values  of  the  vertical  gradients  are 
qualitatively  reasonable.  The  gradients  become  stronger  as  f  decreases,  for  a  given  value  of  a. 

Fig.  9  shows  the  variations  of  the  fluxes  of  potential  temperature  and  moisture  with  height, 
from  the  surface  to  the  PBL  top.  As  before,  the  solid  line  is  for  f  =  1  second,  and  the  dashed  line 
is  for  f  =  10  seconds.  Larger  T  favors  more  linear  flux  profiles. 

As  discussed  earlier,  there  is  reason  to  believe  that  a  is  often  close  to  1/2.  For 
Spi^l  Sp*  «  1,  which  can  be  interpreted  [see  (5.8)  or  Fig.  7]  as  a  close  to  1/2,  (5.12)  and  (5.1 1) 
can  be  approximated  by 


(5.15) 

respectively.  To  obtain  (5.15)  from  (5.11),  we  have  used  (5.8).  In  (5.14)  we  have  kept  terms  of 
first  order  in  Spi^/dp*,  while  in  (5.15)  we  have  kept  terms  of  second  order.  Since  (5.14)  is 
quadratic  in  p ,  differentiating  it  twice  with  respect  to  p  gives  a  constant,  which  is  compatible  with 
our  assumed  constant  value  oi  d\tf  !  dp. 

According  to  (5.15),  dif/  !  dp  vanishes  for  a  =  1/2,  since  this  is  the  limit  dp^  I  dp^^O. 
In  this  same  limit,  (5.14)  reduces  to  the  familiar  linear  flux  profile.  We  conclude  that  our  model 
reduces  to  the  classical  well  mixed  layer  for  <t=  1/2. 


Suppose  that  is  independent  of  height,  as  in  the  water-vapor  example  given  above. 
Then  (5.15)  reduces  to 
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dyf  _  gi 
dp  j 

(5.16) 

For  Fif/^  >  0,  as  in  the  case  of  an  upward  moisture  flux,  (5.16)  predicts  that  ip  decreases  upward. 
This  is  consistent  with  observations  showing  that  the  mixing  ratio  of  water  vapor  typically 
decreases  upward  in  strongly  entraining  boundary  layers  (e.g.,  Wyngaard  and  Brost,  1984).  As  a 

second  special  cas  ',  suppose  that  j(^v/)J«j(^v')5  ’  '^^ich  is  typically  true  for  the  potential 
temperature  flux  in  convective  PBLs.  Then  (5.15)  reduces  to 


dp  a{I  -  a){SpMf[  Sp,  )[ 


a{l-a){Spf4f 


(5.17) 


For  the  case  of  an  upward  surface  potential  temperature  flux,  with  o  <  |  (i.e.,  5pm/6p*  >  0), 


(5.17)  predicts  that  the  potential  temperature  increases  upward,  i.e.  the  flux  is  countergradient; 
again,  this  is  consistent  with  observations  (e.g.,  Wyngaard  and  Brost,  1984). 


If  we  truncate  (5.15)  at  first  order  in  !  dp*,  substitute  into  (5.14),  and  use  (5.9),  we 

obtain 


I  {p-PB){Ps-pW-^)d\(/ 

^  gf  dp' 


(5.18) 


The  first  line  on  the  rhs  of  (5. 18)  is  the  linear  flux  profile  characteristic  of  a  well  mixed  layer,  and 
the  second  represents  a  downgradicnt  diffusion.  According  to  (5.18),  when  the  PBL  fluctuates 
away  from  a  well  mixed  structure,  the  flux  profile  changes  in  such  a  way  as  to  damp  the 


fluctuation.  This  seems  quite  plausible. 
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6.  Summary  and  conclusions 


We  have  presented  what  we  believe  to  be  an  internally  consistent  PBL  model  that  has  a 
parameterized  vertical  structure  but  also  includes  a  simple  second-order  closure. 

In  Section  2,  we  proposed  a  new  way  to  determine  the  updraft  and  downdraft  properties 
applicable  to  a  mass  flux  model  by  using  vertical  motion  statistics  obtained  from  observations, 
LES,  or  even  higher-order  closure.  This  method  cannot  be  used  to  determine  <7  and  Me  in  a  bulk 
mass  flux  model,  however,  because  the  needed  inputs  are  not  available  in  such  a  model. 

In  Section  3,  we  used  imposed  consistency  conditions  among  the  various  fluxes,  and 
applied  suitable  "boundary  conditions"  on  the  updraft  properties  at  the  top  of  the  ventilation  layer 
and  the  downdraft  properties  at  the  base  of  the  entrainment  layer,  to  develop  simple  relationships 
among  mass  fluxes,  fractional  areas,  and  mixing  parameters. 

In  Section  4,  we  presented  methods,  applicable  to  the  bulk  mass  flux  model,  to  determine 
the  convective  mass  flux  and  <y.  So  far  as  we  are  aware,  the  results  obtained  in  Section  4  represent 
the  first  physically  based  method  to  determine  a. 

By  applying  the  mass  flux  model  to  the  interior  of  the  PBL,  we  obtained,  in  Section  5,  an 
idealized  but  physically  based  approach  to  relaxing  the  “well  mixed”  assumption  that  has  been  such 
a  mainstay  of  bulk  boundary  layer  models  over  the  past  25  years  or  so.  In  particular,  the  variance 
budget  equation  determines  the  turbulent  fluxes  and  the  gradients  of  the  mean  state  in  the  interior  of 
the  PBL.  These  vertical  profiles  arc  solved  for  analytically.  Under  appropriate  conditions,  the 
variance  budget  equation  can  reduce  to  downgradient  diffusion,  or  to  a  “compensating  subsidence" 
formula  similar  to  that  used  in  cumulus  parameterizations.  This  represents  a  fir^t  attempt  to  marry 
the  approaches  of  higher-order  closure  and  convective  mass  flux  closure.  Appendix  A  extends  this 
higher-order  closure  approach  to  the  entrainment  and  ventilation  layers. 

In  the  future,  we  hope  to  generalize  the  model  to  allow  <t  and  Me  to  vary  with  height, 
possibly  by  following  the  approach  sketched  in  Appendix  A. 

The  variance  balance  equation  presented  in  Section  5,  together  with  the  TKE-prediction 
equation  given  in  Appendix  B,  comprise  what  amounts  to  a  simple  “second  order  closure,” 
formulated  in  terms  of  the  mass  flux  model,  and  applicable  to  the  bulk  boundary  layer  model.  This 
approach  can  be  put  into  perspective,  as  follows. 
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In  recent  years,  two  new  approaches  have  emerged  for  including  the  effects  of  the  planetary 
boundary  layer  (PBL)  in  large-scale  models.  The  first  involves  coupling  the  large-scale  model  with 
a  "bulk"  PBL  model.  Advantages  of  the  bulk  approach  are  its  simplicity  and  computational 
economy.  A  disadvantage,  up  to  now,  is  its  inability  to  represent  the  internal  structure  of  the  PBL. 

A  second  approach  is  to  make  use  of  a  higher-order-closure  model  in  which  one  or  more 
turbulence  variables  are  prognostically  determined.  This  idea  has  aroused  widespread  interest,  but 
has  been  adopted  in  practice  only  by  K.  Miyakoda's  group  at  the  Geophysical  Fluid  Dynamics 
Laboratory  (Miyakoda  and  Sirutis,  1977;  Miyakoda  etal.,  1983).  Advantages  of  this  approach  are 
its  relatively  high  degree  of  physical  sophistication  and  its  ability  to  predict  the  internal  structure  of 
the  PBL.  Disadvantages  are  its  requirement  for  high  vertical  resolution,  and  its  relative  complexity. 
Both  of  these  lead  to  considerable  computational  expense. 

The  present  model  can  be  viewed  as  an  attempt  to  merge  these  two  approaches,  retaining 
the  advantages  of  each,  and  giving  rise  to  what  we  call  a  "second-order  bulk  model",  or  S.O.B. 
Our  model  also  makes  use  of  a  convective  mass  flux,  a  concept  which  has  been  important  in  the 
development  of  cumulus  parameterizations.  Fig.  10  summarizes  the  conceptual  pedigree  of  our 
model,  relative  to  earlier  models  used  in  boundary-layer  and  cumulus  parameterizations. 

Although  the  approach  presented  here  is  highly  idealized,  it  is  potentially  useful  for 
parameterization  and  also  for  physical  understanding  of  the  results  of  more  complex  models. 
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Appendix  A 


Scalar  Variance  Budgets 
For  the  Ventilation  and  Entrainment  Layers 

In  Section  3,  we  derived  constraints  on  Me  and  <t  by  forcing  the  convective  mass  flux 
model  be  be  consistent  with  the  fluxes  of  y/dne  to  entrainment  and  ventilation.  We  now  derive 

additional  constraints  by  use  of  the  conservation  principle  for  yr'^ ,  as  applied  to  the  ventilation  and 
entrainment  layers. 

Because  ijf  varies  rapidly  across  the  top  of  the  ventilation  layer  and  the  base  of  the 
entrainment  layer,  there  is  vigorous  variance  production  in  these  regions,  but  there  is  also  rapid 
variance  dissipation  by  small  eddies.  The  variance  budget  for  the  ventilation  layer  can  be 
approximated  by 


2p£^pdz  =  2( Vs.  -  Vs)  -  (P^'  ¥Y)s  • 

'*1- 


TA.1) 

The  Ihs  of  (A. I)  represents  dissipation.  The  first  term  on  the  rhs  represents  gradient  production, 
and  the  second  represents  downward  transport  into  the  ventilation  layer  by  triple  correlations.  We 
have  neglected  the  terms  representing  the  local  time  rate  of  change  and  advection  by  the  mean  flow. 


Similarly,  the  variance  budget  for  the  entrainment  layer  can  be  expressed  as 
2  f e^p  dz  --(F^)^Av  +  (PwvY)b  -  V'  ^>8+  2  f '**  v'K  ■ 

(A.2) 

The  Ihs  and  the  first  two  terms  on  the  rhs  of  (A.2)  are  closely  analogous  to  those  of  (A.l).  In  the 
gradient  production  term  of  (A.2),  the  minus  sign  appears  because  of  the  definition  of  Av  -  The 
factor  of  two  that  appears  in  the  production  term  of  (A.l)  is  not  present  in  the  corresponding  term 
of  (A.2),  because  although  the  turbulent  flux  of  ^is  nearly  constant  across  the  ventilation  layer,  it 
drops  from  to  zero  across  the  entrainment  layer,  so  that  its  average  value  for  the  entrainment 
layer  is  half  of  (Fy)B.  The  second  term  of  (A.2)  represents  the  transport  of  variance  into  the 
entrainment  layer  by  triple  correlations.  The  fourth  term  represents  the  rate  at  which  scalar  variance 
is  provided  to  the  newly  entrained  air.  The  fifth  term  represents  variance  production  due  to 
fluctuations  of  5y^in  the  entrainment  layer.  There  is  a  close  analogy  between  (A.2)  and  (3.5). 
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Recall  from  Section  2  that  the  plume-scale  variance  of  y/  is  given  by 


and  that  the  plume-scale  variance  transport  can  be  written  as 


(A.3) 


pw'y/Y=  -2c)i¥u  -  ¥d)^- 

(A.4) 

For  the  LES  results  shown  in  Fig.  3,  a  is  so  close  to  1/2  at  both  levels  S  and  B  that  the  only  safe 
conclusion  is  that  variance  production  and  dissipation  are  closely  balanced  in  both  the  entrainment 

and  ventilation  layers.  If  =  [pw'\ff'\j/')g>  the  variance  flux  passes  through  the  interior 

of  the  PBL  without  convergence  or  divergence. 


These  ideas  may  have  some  relevance  to  mesoscale  cellular  convection,  if  our  mass  flux 
model  can  be  applied  to  the  mesoscale  convective  circulations.  For  the  case  of  open  cellular 
convection,  surface  heating  is  believed  to  dominate  (e.g.,  Arakawa,  1975)  and  observations 
suggest  that  as  «  /.  According  to  our  model,  the  entrainment  layer  with  ctb  «  I  is  dissipating 
more  variance  than  it  produces.  This  seems  physically  plausible  for  open  cells.  For  closed  cellular 
convection,  cloud-top  cooling  is  believed  to  be  dominant  (Arakawa,  1975)  and  observations 
suggest  that  cTs  =  /  .  According  to  our  mottel,  the  entrainment  layer  with  cr^  =  /  is  producing  more 
variance  than  it  is  dissipating.  This  seems  physically  plausible  for  closed  cells. 


We  now  define  nondimensional  measures  of  the  dissipation  rates,  denoted  by  kv  and  ks  . 
respectively; 


3  kyV a^d  -  (Ts)( ¥.  -  ¥<t)l  - 

2  r**£^  dz  3  k^E  CgU  -  (7g)(  -  v^j)*; . 


(A.5) 

(A.6) 


Rapid  variance  dissipation  rates  in  the  ventilation  and  entrainment  layers  correspond  to  large  values 

of  ky  and  ks ,  respectively.  Note  from  (A.3)  that  the  rhs  of  (A.5-6)  are  proportional  to  \|/^  for  the 
ventilation  and  entrainment  layers,  respectively.  The  definitions  of  ky  and  ks  ,  given  by  (A.5-6). 
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are  motivated  by  the  idea  that  the  rates  of  variance  dissipation  in  the  ventilation  and  entrainment 
layers  should  be  related  to  the  actual  values  of  the  variances  at  the  edges  of  those  layers. 

Using  the  results  of  Section  3  with  (A.4-5),  we  rewrite  (A.l)  as 


<T5  = 


2-Xv 

kvXv-^^U-Xv)  ' 


(A.l) 

All  referen''T  to  yrhas  dropped  out  of  (A.7).  indicating  that  kv  is  independent  of  species. 

Now  turning  to  the  entrainment  layer,  we  can  use  (A. 3),  (A. 4),  and  (A.6)  with  the  results 
of  Section  3  to  rewrite  (A. 2)  as 


( ¥u  -  ¥d)BXE^c.B{-W  -  <^b)  +Xe\xe^b  +  -(1-  era)} 

Ilf  „  _  _ 

+( ¥u  -  ¥'S'^  dz  =  0. 

,  (A.8) 

We  require  that  (A.8)  hold  for  all  species,  including  those  for  which  Sy/is  zero.  It  follows  that 


Cb  = 


IceXe-Xe-^Xe 

I-^^eXe-^Xe-^Xe 


(A.9) 


Again,  ail  reference  to  v^has  dropped  out,  showing  that  kg  is  independent  of  species. 

Since  (A.8  -  9)  must  apply  for  all  species,  including  those  for  which  Sy^is  not  zero,  we 
conclude  that 


Jz.  E  i)Jz.  ^ 

(A.  10) 

Using  (A.  10),  we  can  determine  the  rate  at  which  fluemations  of  in  the  entrainment  layer 
generate  fluctuations  of  v' there. 

Using  the  results  derived  above,  we  can  show  that  the  ratio  of  dissipation  to  production, 
for  the  ventilation  layer,  satisfies 
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(A.ll) 


y^s-~  .  2(1  —  0^) 


The  corresponding  ratio  for  the  entrainment  layer  is 


~L  ^  ,  r  Xe{1-^^b) 

2  f  IjFsT  dz  i(^B-XE(^-<yB).' 

Jzm  ^ 


(A.  12) 

Note  from  (A.  1 1)  and  (A.  12)  that,  since  Xv  and  Xe  are  both  on  the  order  of  lO'^,  and  if  <7  is  close 
to  1/2,  dissipation  and  production  must  nearly  balance  in  the  ventilation  and  entrainment  layers. 
Such  near  balances  arc  observed  (e.g.,  Caughey  and  Palmer,  1979),  and  have  been  predicted 
through  LES  (see  for  example  Fig.  7  of  Mocng  and  Wyngaard,  1989). 

Now  combine  (A.7)  with  (3.4),  to  obtain  a  quadratic  equation  for  Xv  ■  The  solutions, 
which  arc  both  physically  relevant,  arc  given  by 

2G,  ^kv±  ^14  -8Giky  +  4G} 


2kvGv 


Using  (3.4),  we  can  obtain  the  corresponding  solutions  for  <75 : 


(A.13) 


_  {^ky  —  2Cy)  T  ■yj ky  ~  8Gyky  +  4Gy 

Uy  ■ 

(A.  14) 

The  condition  that  the  discriminant  in  (A.13)  and  (A.  14)  be  non-negative  implies  that  /:v/must  be 
sufficiently  large,  for  a  given  Gv  ■  The  discriminant  vanishes  for 


ky  =  {ky)^  =  4Gi^  +  2Gy^4Gf-l . 


(A.  15) 


For  convenience,  let 


(/tv)  ■ 

'  ^  'mm 


(A.  16) 
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Fig.  1 1  shows  plots  of  Xw  and  (Js  against  Gy  and  ry  .  Both  solutions  are  shown.  Note  that,  for 
sufficiently  large  Gy,  <Js  depends  only  on  rv.  For  Gy  >>  I,  the  following  approximate  relations 
hold: 


From  (A.  18),  we  conclude  that 


=  8Gy  , 


<Ts  = 


Xv  = 


(A.  17) 

fy  T  -Jrj/lryr  —  I) 

2ry 

(A.18) 

ry±^ry{ry-l) 

2ryGy 

(A.  19) 

-  ^  tTc  ^  /  with  the  choice  of  the  minus,  or 
2  ^ 


0  ^  cr5  S  -j  with  the  choice  of  the  plus. 


(A.20) 


By  substituting  (A.  17  -  19)  back  into  (A.l)  and  (A.5),  we  can  show  that  the  approximation 
Gv>>l,  on  which  (A.17-19)  are  based,  corresponds  to  a  balance  between  dissipation  and 
production. 


A  parallel  analysis  for  the  entrainment  layer  leads  to 


Xe  = 


2GE  +  kE  +  l±  kl^]ki  +  {2-  4Gl)kE  1 
2G£{1  +  k^) 


{A.21) 


and 


Ob  = 


2GE  +  kE  + It  k^^Jkj  +  (2-  4Gl)k£  + 1 

2(1 +  kE) 


(A.22) 

Again,  both  solutions  are  physically  relevant  The  condition  that  the  discriminant  in  (A.21  -  22)  be 
non-negative  implies  that  k£  must  be  sufficiently  large.  The  discriminant  vamshes  for 
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(A.23) 


Let 


_  ^E 
V  *-/min 


(A.24) 

Fig.  13  shows  plots  of  Xe  against  Ge  and  r£.  The  situation  is  very  siniilar  to  that  found  for 

the  ventilation  layer.  For  sufficiently  large  Ge,  Ob  appears  to  depend  only  on  rg  •  We  can 
demonstrate  that,  for  G^  >>  /,  the  following  approximate  relations  hold: 


<Tg  = - - - 

2^£ 


(A.25) 


(A.26) 


_  r£±V'£('£~^) 

— w— 

(A.27) 

Note  that  (A.26  -  27)  closely  parallel  (A.  19).  Fig.  12  shows  Ob  as  a  function  of  rg ,  and/or  as  as  a 
function  of  ry,  as  given  by  (A.  18)  and/or  (A.26). 

From  (A.26),  we  can  conclude  that 


—  ^Ob^I  with  the  choice  of  the  plus,  and 
2 


O^Cb^-^  with  the  choice  of  the  minus. 


(A.28) 

The  main  results  of  this  Appendix  are  (A.  18- 19)  and  (A.26- 27).  Our  expressions  for  cr5 
and  Ob  automatically  satisfy  0  <  <t  <  /.  If  ry  and  rg  were  known,  we  could  find  05,  Os,  and  Xe 
and  Xe-  The  values  of  Aff.5  and  McJb  would  follow  immediately  from  (3.4)  and  (3.9),  respectively. 
The  interior  value  of  Me  could  be  obtained  by  an  approach  analogous  to  that  used  in  Section  4.  We 
would  thus  be  in  a  position  to  allow  a  to  vary  linearly  with  height,  and  Me  to  vary  quadratically. 
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This  would  allow  improved  agreement  with  the  LES  results  presented  in  Fig.  3.  Pursuit  of  this 
generalization  of  our  model  is  left  to  the  future. 

APPENDIX  B 

A  Method  to  Predict  the 

Vertically  Integrated  Turbulence  Kinetic  Energy 

The  vertically  integrated  conservation  law  for  the  TKE  can  be  written  as 

g'Sf>„^*Ee„  =  B  +  S-D. 

(B.l) 

Here  B,  S.  and  D  represent  production  by  buoyancy,  production  by  shear,  and  dissipation, 
respectively,  and  dpi^  is  the  pressure-thickness  of  the  PBL. 

The  vertically  integrated  dissipation  rate  and  the  vertically  averaged  TKE  are  assumed  to  be 
related  by 

D=p„(e„la,f\ 

(B.2) 

where  py,  is  the  vertically  averaged  PBL  density,  and  a,  =0.163  is  a  dimensionless  constant. 

The  buoyancy  production  integral,  B,  is  of  the  form 

B^Kr^dp. 

p 

(B.3) 

where  k  is  Poisson's  constant  This  integral  can  be  evaluated  using  the  methods  of  Randall 
(1984). 

The  shear  production  integral,  5,  is  of  the  form 
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(B.4) 


9V 

dp 


dp. 


We  divide  this  integral  into  three  parts:  shear  production  in  the  surface  layer,  in  the  interior  of  the 
PBL,  and  in  the  entrainment  layer.  The  surface  layer  shear  production  is  approximately  given  by 

C 

(B.5) 

Here  the  second  equality  follows  from  the  assumption  that  the  surface  stress  in  parallel  to  the 
surface  wind. 

The  interior  shear  production  can  be  evaluated  by  straightforward  methods  once  the  profiles 
of  the  stress  and  the  wind  are  known. 

Finally,  the  entrainment  layer  shear  production  is  approximately  given  by 

syy  •  4- = iFy)i\Av\ = -tcK. 

(B.6) 

In  (B.6),  the  factor  of  1/2  arises  because  Fv  decreases  from  {Fv)b  to  zero  across  the  entrainment 
layer,  so  that  its  average  value  inside  the  entrainment  layer  is  {Fv)b  I  2.  The  second  equality 
comes  from  applying  (3.5)  to  the  momentum  budget  of  the  entrainment  layer. 

APPENDIX  C 

Motivation  for  (5.8  -  9) 

Dropping  the  time-dependence  of  (5.3)  and  simplifying,  we  obtain 

(/-2<y)<^^y  <T(l-a)  dip 

8'^disl^c 

(C.l) 
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This  is  basically  the  same  as  (5.7),  except  that  we  have  not  introduced  dp*.  We  solve  (C.  1)  as  a 
first-order  ordinary  differential  equation  for  Fy, ,  assuming  constant  coefficients  and  a  constant 
value  of  d\\f  !  dp.  The  solution  is 


P  ^  ^  r  an-  ajp  -Ps) 

^  o(l~a)dp 


(C.2) 

where  A  is  the  constant  of  integration.  Notice  that  the  quantity  (I  -  2(7)  appears  in  the  denominator 
of  the  argument  of  the  exponential.  As  <T  passes  through  1/2,  the  argument  of  the  exponential 
diverges  to  plus  and  minus  infinity  —  unacceptable  behavior.  We  cannot  solve  this  problem  by 
choosing  A  =  0,  because  then  Fy,  would  be  independent  of  height  and  could  not  satisfy  its 
boundary  conditions  at  the  surface  and  the  PBL  top. 

Our  interpretation  is  that  must  depend  on  (7.  This  is  acceptable,  since  the  a- 
dependence  of  is  not  known  a  priori.  We  choose  the  (7-  dependence  of  in  such  a  way 
that  Fyf  varies  continuously  with  <7. 

This  we  can  do  by  taking  proportional  to  (/  -  2(j),  i.e. 


i  =  {l-2(7)Xji,, 


where  f  is  a  constant  (i.e.,  independent  of  o).  Then  (C.2)  is  replaced  by 


(C.3  a) 


(C.4a) 


where 


dp*  = 

o(l-<7) 


(C.5  a) 

Inspection  of  (C.4a)  shows  that  the  flux  does  indeed  vary  continuously  with  cr.  This  approach  is 
unacceptable,  however,  because  (C.3  a)  will  inevitably  give  negative  and  hence  physically 
impossible  values  of  for  some  values  of  <7. 


We  can  ensure  non-negative  and  also  make  vary  continuously  with  a  by  choosing 


37 


x^{l-2afx^„. 


(C.3  b) 


with  X  positive  and  independent  of  a.  The  solution  of  (C.l)  is  then 


(C.4b) 


(C.5  b) 

This  is  the  line  of  reasoning  that  led  us  to  adopt  (5.9),  with  r  independent  of  a. 

An  implication  of  (5.9)  is  that  becomes  infinite,  i.e.  the  variance  dissipation  rate  goes 
to  zero,  for  <t  =  1/2.  Of  course,  this  is  at  best  an  idealization  of  how  nature  works.  From  the  LES 
results,  we  estimate  that  (/  -  2<t)  ~  10**  or  lO*^,  which,  together  with  the  results  shown  in  Fig.  6. 
suggests  that  for  this  particular  case  f  is  on  the  of  order  one  to  ten  seconds. 

Why  should  x^j  increase  as  a  approaches  1/2?  A  possible  physical  interpretation  is  that 
the  largest  (absolute)  convective  vertical  velocities  increase  as  a  departs  from  1/2.  For  example,  if 
a  «  1,  the  updraft  vertical  velocities  must  be  relatively  vigoious,  implying  strong  lateral  shears 
and  suggesting  strong  lateral  mixing  and  variance  dissipation.  Similarly,  for  o  close  to  one  the 
downdraft  vertical  velocities  must  be  relatively  vigorous.  The  intensity  of  the  individual  drafts  is 
minimized  for  a  =  1/2,  suggesting  that  the  dissipation  rate  is  minimized  then  too. 

We  must  also  ask  why  T  is  so  shon.  From  (5.2)  and  (5.9),  the  dissipation  rate  can  be 
written  as 


=  Mc\ 


^  ^  .  ...J P-Pb^ 


{I -2a)  dp 


+  A  exp 


6fK 


where 


5fu  s 


gM,x 


a{I-a){l-2a) 


We  have  used  (C.3  b)  and  (C.5  b)  in  Section  (5). 


£^={l-2a)‘ 


V  ^  / 


(C.6) 
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We  can  interpret  f  as  the  intense  local  dissipation  rate  in  those  ponions  of  the  convective 
circulations  where  dissipation  actually  occurs,  and  {I-  2 erf  as  the  fraction  of  the  domain  within 
which  dissipadon  is  localized.  The  volume-averaged  dissipadon  rate  is  then  given  by  (C.6). 
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Table  1  Summary  of  observations  and  simulations  of  the  fractional  area  covered  by  rising 
motion  in  several  studies.  Here  {  )  denotes  the  expected  value. 


Author 

Source  of  Data 

Sampling  Method 

Used  to  Define 
“Updraft” 

Mean  o 

Manton  (1977) 

Observations 

T  >  Tth  such  that 

0.42 

Coulman  (1978) 

Observations 

Similar  to  Manton  (1977) 

0.38 

Lamb  (1978) 

LES  results  for  clear 
PEL 

Vertical  velocity  is  positive 

0.45 

Greenhut  and 
Khalsa(1982) 

Observations 

Vertical  velocity  exceeds  a 
positive  threshold.  Three 
domains  are  defined;  One 
for  updrafts,  one  for 
downdrafts,  and  a  third 
“environment.” 

0.16 

Lenschow  and 
Stephens  (1982) 

Observations 

Moisture  fluctuation  is 
positive 

0.25 

Young  (1988  a) 

Observations 

Vertical  velocity  is  positive 

0.46 

NichoUs  (1989) 

Observations  of 
stratus-topped  PEL 
over  the  North  Sea 

Vertical  velocity  is  positive 

0.70 

Schumann  and 
Mocng  (1991  a) 

LES  results 

Vertical  velocity  exceeds  a 
positive  threshold.  Three 
domains  are  defined:  One 
for  updrafts,  one  for 
downdrafts,  and  a  third 
“environment” 

0.35 
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Table  2 


The  LES  entrainment  mass  flux  in  lO’^  kg  m'^  s'^  as  determined  by  methods  HA. 
HB,  QA,  QB,  and  Z.  See  text  for  explanation  of  the  methods. 


History  Record  Number 


1 

2 

3 

4 

5 

6 

Average 

HA 

0.803 

0.768 

0.816 

0.814 

0.797 

0.841 

0.806 

HB 

0.660 

0.639 

0.666 

0.691 

0.741 

0.767 

0.694 

QA 

0.819 

0.786 

0.859 

0.835 

0.853 

0.833 

0.831 

QB 

Z 

0.633 

0.625 

0.670 

0.671 

0.760 

0.733 

0.682 

0.656 

Table  3  Values  of  ;|^v  and  ;|<£  determined  from  the  LES  results  by  various  methods. 

Subscript  q  indicates  that  we  have  used  y/  =  qt,  while  subscript  h  indicates  that  we 
have  used  y/-  h.lti  the  case  of  x'v,h  we  have  taken  into  account  the  weak  radiative 
forcing  of  the  ventilation  layer.  This  slightly  reduces  the  values  of  Xv.h  .  bringing 
them  into  better  overall  agreement  with  those  of  xy.q  • 


History 

Record 

Number 

1 

2 

3 

4 

5 

6 

Average 

0.0124 

0.0124 

0.0130 

0.0121 

0.0125 

0.0119 

0.0124 

Xv.h 

0.0126 

0.0137 

0.0140 

0.0130 

0.0140 

0.0127 

0.0133 

XE..J 

0.0097 

0.0104 

0.0101 

0.0094 

0.0109 

0.0106 

0.0101 

Xz.h 

O.OIOI 

0.0106 

0.0101 

0.0097 

0.0107 

0.0110 

0.0104 

X’v,  h 

0.0113 

0.0124 

0.0128 

0.0117 

0.0127 

0.0114 

0.0121 
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FIGURE  LEGENDS 


Figiire  1:  Diagram  illustrating  the  assumed  structure  of  the  PBL.  The  interior,  which  is 
represented  by  two  layers,  is  bounded  above  by  a  thin  entrainment  layer  and  below 
by  a  thin  ventilation  layer.  Convective  circulations  occur,  with  rising  branches 
occupying  fractional  area  <J.  The  ascending  and  descending  branches  have  different 
thermodynamic  soundings  and,  therefore,  different  cloud  base  levels. 

Figure  2:  The  vertical  profiles  of  h  and  q, ,  plotted  as  functions  of  height. 

Figure  3:  a)  The  vertical  profiles  of  and  obtained  by  the  method  described  in  Section  2 
(solid  curves),  and  also  by  direcdy  sampling  the  LES  results  (dashed  curves),  b)  The 
vertical  profiles  of  a  (solid  curve)  and  Me  (dashed  curve),  plotted  as  functions  of 
height  normalized  by  the  PBL  depth.  Values  are  not  plotted  above  the  PBL  top,  since 
trand  Me  have  little  meaning  there. 

Figure  4:  Profiles  of  \ffu  and  calculated  directly  by  conditional  sampling  (based  on  the  sign 
of  the  vertical  velocity)  and  from  the  fluxes  using  (2.10  - 1 1).  a)  for  h  /  Cp\  and  b)  for 
qt .  The  solid  lines  show  the  updraft  and  downdraft  properties  obtained  from  (2.10  - 
11),  and  the  dashed  lines  show  the  corresponding  values  based  on  conditional 
sampling. 

Figure  5:  For  the  LES  results,  the  ratio  of  the  subgrid-scalc  flux  moisture  flux  to  the  total  flux. 
The  ventilation  and  entrainment  layers  are  indicated  by  stippling. 

Figure  6  Vertical  profiles  of  as  determined  from  the  LES  results  for  dissipation  of  h  (solid 
line)  and  qt  (dashed  line). 

Figure  7  The  variation  of  Spi^  /  Sp*  with  cr,  for  the  particular  case  Me  =  0.4  kg  m‘2  s‘^  Sp^  = 
50  mb.  The  solid  line  is  for  f  =  1  second,  and  the  dashed  line  is  for  r=  10  seconds. 

Figure  8  As  in  Fig.  7.  but  for  the  variation  of  dip  I  dp  with  a:  a)  potential  temperature,  in 
units  of  K  (50  mb)-^;  b)  the  mixing  ratio  of  water  vapor,  in  units  of  g  kg-^  (50  mb)-f 
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Figure  9 


Figure  10 


Figure  1 1 


Figure  12 


Figure  13 


As  in  Fig.  7,  but  for  the  variations  of  the  fluxes  of  potential  temperature  and  moisture 
with  height,  from  the  surface  to  the  PBL  top:  a)  Potential  temperature  flux  for  a  = 
0.52;  b)  moisture  flux  for  a  =  0.52;  c)  potential  temperature  flux  for  a  =  0.48;  d) 
moisture  flux  for  G  =  0.48. 

Diagram  summarizing  the  relationship  of  the  present  model  tc  earlier  models  used  in 
boundary-layer  and  cumulus  parameterizations. 

Plots  of  Xv  and  Os  against  ky  and  ry,  as  given  by  (4. 13)  and  (4. 15).  Panels  a  and  b 
for  for  the  choice  of  the  plus  sign  on  the  discriminant,  and  panels  c  and  d  are  for  the 
choice  of  the  minus  sign. 

Plots  of  Xe  and  Gb  against  Jtf  and  r^,  as  given  by  (4.21  -  22).  Panels  a  and  b  for  for 
the  choice  of  the  plus  sign  on  the  discriminant,  and  panels  c  and  d  are  for  the  choice 
of  the  minus  sign. 

Plot  of  (Tg  as  a  function  of  r£ ,  and/or  of  Os  as  a  function  of  ry. 
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